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Abstract: The present paper analyze the constraints on the confluent Heun type-equation, (aspr^-l- 
a 3 , 2 r)y'' + (a 2 ,or^ +a 2 ,ir + a 2 , 2 )y' -(Tifir + Ti^i)y = 0, where |a 3 ,ip + |a3,2|^ * 0, and aij,i = 3,2, l,j = 
0,1,2 are real parameters, to admit polynomial solutions. The necessary and sufficient conditions 
for the existence of these polynomials are given. A three-term recurrence relation is provided to 
generate the polynomial solutions explicitly. We, then, prove that these polynomial solutions are 
a source of finite sequences of orthogonal polynomials. Several properties, such as the recurrence 
relation, Christoffel-Darboux formulas and the moments of the weight function, are discussed. We 
also show a factorization property of these orthogonal polynomials that allow for the construction 
of other sequences of orthogonal polynomials. For illustration, we examines the quasi- exactly 
solvability of the (p, g)-hyperbolic potential V(r) = -Vo sinh^(r)/cosh'^(r), Vo > 0,p >0,q>p. The 
associated orthogonal polynomials generated by the solutions of the Schrodinger equation with the 
(4,6)-hyperbolic potential are constructed. 
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I. INTRODUCTION 


The second-order differential equation 


(Py 


dy 


(o3,ir + 03, 2 r) — + (02,oc + 02,ir -f 02 , 2 ) — - ('ri,or + Ti^i)y = 0 


dr^ 


dr 


®3,i ®3,2 


®2,o U,o 


( 1 ) 


with parameters Qij € R, i = 3,2, 1, j = 0,1,2 has two regular singular points r = 0 and r = -03,2/03,1 with exponents 
{0,1 - 02,2/03,2} and {0,1 - 02,1/03,1 + 02,2/03,2 + 02,003,2/0! i), respectively, in addition to an irregular singular point 
at infinity unless 02,2 = U.! = 0. 

Although equation ([T]) generalized the standard confluent Heun equation [l|, we shall keep the same name ‘the 
confluent Heun equation’ for Jp as well. The name confluent driven from the coalescence process of singularities of 
the general Heun equation [ll-Q. 

Over the past few decades, the confluent Heun equation appeared intensively in mathematical physics applications. 
Including but not limited to, in studying the spherical Coulomb functions, the spheroidal wave functions, the Mathieu 
equation and confinement potentials in cmantum mechanics. Most recently, it also found some interesting applications 
in General Relativity and Astrophysics |4l4l0l|. 

The purpose of this article is twofold. First, we analyze the necessary and sufficient conditions under which the 
differential equation m admit the polynomial solutions yn{r) = T,k=o^kP, n = 0,1,2 ,..., Where we set up a scheme 
for constructing the polynomial coefficients Ck explicitly. Second, we establish a correspondence between each solution 
?/„(r) and a finite sequence of orthogonal polynomials {Tfe(Ti,i)}^^Q. We discuss some properties of these orthogonal 
polynomials such as the recurrence relation, factorization, Christoffel-Darboux formulas, and the moments of the 
weight function. 

As an illustration of these polynomials, we explore the quasi- exact solutions of the one-dimensional Schrodinger 
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equation with a hyperbolic potential 


h'^Cp sinh^(a:/d) r \ 

2 m dx^ cosn^(x/a) 


oo < x < oo, 


( 2 ) 


with special attention to the case of p = 4 and q = 6 . We show under certain transformations, this equation reduces 
to a confluent Heun equation of type ©■ We discuss the associated sequence of orthogonal polynomials {(Pk}k=o 
generated by the solution ifjmny 1 and investigate some of their properties. 


The article organized as follows. In Section II, we gave the necessary and sufficient conditions for the polynomial 
solvability constraints of equation ([T|) with a general procedure for formulating the polynomial solution explicitly. In 
Section III, we establish a correspondence between that the solution ?/„ and a set of orthogonal polynomials {^k}k=o- 
In that section, we also derive the recurrence relation, the Christoffel-Darboux formula and prove the factorization 
property of these polynomials. The weight functions and some explicit expressions calculating the moments of the 
weight functions also discussed. In Section IV, we discuss the quasi- exact solvable quantum model, given by a 
hyperbolic potential and explicitly construction the finite sequences of orthogonal polynomials associated with the 
wavefunction solutions. 


II. POLYNOMIAL SOLVABILITY CONSTRAINTS 


The n-th derivatives of equation m yields 

(03,+ 03,+ (02,0?’^ + [02.1 + 27103,i]r + ^^03,2 + 02,2)2/*'"’^^^ 

+ ((2no2,o - Ti,o)r + n{n - 1)03,1 + na2,i - ti,i) + (n(n - 1)02,0 - oryo) = 0 , n = 0 , 1 ,- ( 3 ) 

Thus, the necessary condition for the n-degree polynomial solution is ti,o = na2fi,n = 0 , 1 ,... while the sufficient 
condition follow using the recurrence relation 

(fc-H 1 ) (fco3,2 + 02,2) Cfc+i + (fc(fc - 1)03,1 + fco2,i - Ti,i))Cfe + ((fc - 1)02,0 -'ri,o))Cfc_i = 0 , C-i = 0 , Co = 1 , ( 4 ) 

according to the following simple procedure. 

For the zero-degree polynomial solution yo{r) = 1 : the necessary and sufficient conditions, respectively, reads 

n.o = 0 , Ti,i = 0 . ( 5 ) 

For the first-degree polynomial solution 

yi{r) = l+—r, (6) 

02,2 

the necessary and sufficient conditions, respectively, reads 


Ti,o - 02,0, 

For the second-order polynomial solution 


-Tl,l 02,2 

-02,0 02,1 -Ti,i 


= 0 . 


, , Til 202,002,2 - 02.lTl,l H-„ 

y2{r) = l+—r+ --- - -;- —^r , 

02,2 202,2(02,2 + 03,2) 


the necessary and sufficient conditions, respectively, reads 

Tl,0 = 2 o2,o, 


-Ti,i 02,2 0 

-2o2,0 02,1 -Ti,i 2o3,2 + 2o2,2 

0 “02,0 2o3,i -H 2o 2,1 - Ti,i 


= 0 . 


( 7 ) 


( 8 ) 


( 9 ) 


For the third-order polynomial solution 


ysir) = 1 + —^r + 
02,2 


Ti,i 3 o 2,0O2,2 - a2,lTi,i H-rf 2 


202,2(02,2 + 03,2) 


-602,002,2(02,1 + 03 1) + 02 .o( 702,2 + 4o3 2)+! 1 + (02,1 - Ti.i)(2(o 2,1 + O3 i) - Ti i)ti.i 3 
H-^^^ T 

602,2(02,2 + 03,2)( o 2,2 + 203,2) 


( 10 ) 
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the necessary and sufficient conditions, respectively, reads 


Tlfi - 302 , 0 ; 


-''"1,1 02,2 
-302,0 02,1 -Ti,i 
0 —202,0 

0 0 


0 

203,2 + 2 o2,2 

203,1 + 2 o 2,1 — Til 
- 02,0 


0 

0 

603,2 + 302,2 
603,1 + 302,1 — Ti,i 


= 0 . 


( 11 ) 


In general, for an n-degree polynomial solution of the confluent Heun equation, we have the following theorem. 

Theorem II. 1 . The necessary and sufficient conditions for an n-degree polynomial solution of equation o , respec¬ 
tively, are 


Ti,o - riO2,0; n - 0,1,..., 


and the vanishes of the (n+ 1) x (n+ 1) determinant 


( 12 ) 


-Tl,l 

02,2 

0 

0 

0 

0 

-Tl,0 

02,1 - n,i 

203,2 + 2o2,2 

0 

0 

0 

0 

02,0 -U,0 

203,1 + 2o 2,1 — Til 603,2 + 302,2 

0 

0 

0 

0 

202,0 - Tl,0 

603,1 + 302,1 - Ti,i . . 

0 

0 

0 

0 

0 

0 

(n- 1)02,0 - Ti,o 

n{n - 1)03,1 + no2,i - ti,i 


Theorem III. II generalized the two conditions on the standard confluent Heun equation to allow for the n-degree 
polynomial solutions i®- 

Remark II.1: If 03,2 = 02,2 = 0 , the differential equation ([T]) reduces to 




03,H + »'(02,0?' + 02,1)^ - {Tipr + Ti i)y = 0, 
ar^ dr 


that have the n-degree polynomial solutions expressed in terms of the confluent hypergeometric function as 

yn-k{r) = r^iFi ik-n-, 2 k + ^ fc = 0,1,.. .n, 

V 03,1 03,1 ) 

provided that: ti,o = na2,o and ri,i = fc(A;-l)a3,i + /ca2,i, fc = 0,I,2,...,n, n = 0,I,2,_ 


( 14 ) 


( 15 ) 


III. FINITE SEQUENCES OF ORTHOGONAL POLYNOMIALS 

Subject to the necessary and sufficient conditions (USD, the polynomial solutions ([ 5 ]), (|S]) and (HHD of the differential 
equation o can be generalized as 


lir) = E 




k =0 k\ Og 2 


/ 02,2 \ 
\ 03,2 )j, 


n = 0,I,2, 


( 16 ) 


where (a)„ refer to the pochhammer symbol {a)n = a{a+ I)... (a-n + I) satisfies {-n)k = 0 for any positive integers 
k> n. 


The coefficients k = 0 ,... ,n are a sequence of /c-degree polynomials in the parameter ( = ti,i. This sequence of 

polynomials will be our focus in the present work. Indeed, for each n, the solution ynif) of the differential equation 
o generate a finite sequence of polynomials {2fc(Ti,i)}fc=o satisfying a number of very interesting properties. 


A. Three-term recurrence relation 


The first property is that they satisfy a three-term recurrence relation 
^k+iiO + (fc(fc- 1)03,1 H- ka2,i - C) Tjf(C) + fc((fc- 1)03,2 + 02,2) ((fc- 1)02,0 -Ti,o) Tjf_i{C) = 0, 0<fc<n-I, ( 17 ) 
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initiated with 2 ’"i(C) = Oi ^o*(C) = 1 - This recurrence relation can be obtained by the fact that y„(r) satisfies the 
differential equation O- Using the necessary condition ti q = na2.o, equation reduce to 

^k+iiO = 1)03.1 -fca 2 ,i) 2 ’fc(C) -k{k-n-l)a 2 ,o{{k-l)a 3^2 + a 2 , 2 )^k-i{ 0 , ^-liO = <(0 = 1, ( 18 ) 

for fc = 0 , 1 , 2 ,... n - 1 . Here, n is a (fixed) nonnegative integer refer to the degree of the polynomial solution (TlBl) . 

From the general theory of orthogonal polynomials and according to Favard’s Theorem [Iffl), form a 

sequence of orthogonal polynomials. The first few expressions of these polynomials reads 


^liO = c, 

®?(C) = “ <3^2,iC + ^102,0 02,2) 

“ ( 302,1 + 2a3,i)C^ + {^0,2 I + 202,103,1 + 02,002,20 - 2(02,2 + O3,2)(o2,o “ 02,00)) C “ 2no2,002,2(02,1 + 03,1) 

^r(C) = “ 2 ( 302,1 + 4 o 3 ,i)C^ + {lla 2 i + 2602,103,1 + 1203 3 + 202,0(02,2(80 - 4 ) + 03,2(471 - 7 )))((^ 

- 2 (3(02,1 + 03,1) (02,2 - 202,002,2 + 202,103,1) - 02,003,2(802,1 + 603,1) + 002,0(02,2(702,1 + 10o3,i) + 6(02,1 + 03,1)03,2)) C 
+ 37102,002,2(2(02,1 + 03 ,i)(o 2 ,l + 203,1) + 02,0(0 - 2 )(o2,2 + 203,2)). 


An important consequence of the three-term recurrence formula ( 1181 ) is the following Christoffel-Darboux formula. 


Theorem III.l. For ( + C,', 

^ ip;(C)g;(C') gfe",,(C)g,"(CO-gfe"(C)g-fe".i(C') 

.=0 j \ (02,003,2)^- (If )^. (-0), k \ (02,003,2)^= (If (-o)fc(C - CO ’ 

while for C^ ^ C 

^ (gr(c))' [g-fcOi(c)]'gfc(c)-[g-fc(c)]'gfcOi(c) 

.=0 j \ (02,003,2)^ (If )^. (-0), ■ k \ (02,003,2)^= (If (-O)fe 

Proof. The recurrence relations (IT^ for C and reads, respectively, 


fc = 1 , . . . , 71 - 1 , 


k = 1 , . . . , 77 - 1 . 


( 19 ) 


( 20 ) 


2 feVl(C) = C^kiO - [fc(fc- 1)03,1 + fco2,l]2’fc (C) - fc(fc-n - l)o2,o((fc - 1)03,2+ 02,2)2’fc_i(C), 

^k+liC) = C'PkiC) - [fc(fc- 1)03.1+ fca 2 ,l] 2 ’fc(C 0 -fc(fc-0 - l)o2,o((fc - 1)03,2+ 02,2)!Pfe-l(C')- 

Multiplying the first by 2 ’^(CO ^ad the second by 2 ’^(C) then subtruct, the resulting equation reads 

(C-COa’fc (COa’fc (C) = Qfc 7 i(C,C 0 - Afc.igfe(C,CO 

where Afe+i = fc (fc - n - l)o2,o((fc - 1)03,2 + 02,2), and Qk+iiC,C) = ^k+iiO^kiC) - (C) 2 ’fc+i(C 0 - Thus, recursively 

over k we have 


(C - COs’fe (Os’fe (CO = QkriiC, CO - XkriQkiC, CO 
(C - COs'fc-i(C)2’fe-i(CO = QkiC, CO - AfeQfe-i(C,CO 


(C-CO<(C)2’o”(CO = Qi(C,CO- 

From which, it is straightforward to obtain 


(C-CO 


^kiO^kiC) + Afe^i2’fc”_i(C)^Pfe-i(C0 + Afc+iAfe2’fe”_2(C)2’fe-2(C0 + Afc^iAfeAfc_ifffcU 3 (C)C 3 (C 0 + 


+ Afe+iAfcA/c-lA/c-2 ■ . ■ A 2 !P, 5 *(C) 2 o*(C 0 


= Qfc+i(C,CO 


Divide both sides by (C - C0 Afe+iAfcAfe_iAfe_2 ... A2 and sum over k we finally gets 

1-1 2’fcOl(C)2’fe"(CO - iPm^krliC) 




j=o Aj+iAjAj_i... A; 


- (Afc+iAfeAfc_i... A2) 


C-C' 
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Formula (El) then follows using 


fc+1 

Afe+iAfcAfe-i... A2 = n Aj 

J=2 


k\ (02,003,2)^ 





k = 1 ,... ,n, 


and finally the formula (IT^ follows by evaluating the limit of both sides of (flTll as C ^ C- ^ 

Theorem III. 2 . The zeros {Ci}J=o ^ke orthogonal polynomial = 0 ,...,n- 1 are eigenvalues of the 

following tridiagonal matrix 


'’oo 1 ' 

/ 3 i Oil 1 
/32 0-2 1 

f^k-l Oik-1 1 
Pk Ctkl 


( 21 ) 


where aj = j(j - 1)03,1 + ja2,i, for j > 0 and (ij = j {j - n - l)a2,o((j - 1)03.2 + 02,2), for j > 1. 
Proof. We first write the recurrence relation dT 51 ) as 


C 2 >;(c) = i(c) + i(C), J > 0 , 


We now take j = 0 , 1 , 2 ,... fc to form a system with a matrix form 


C^oiO ] 


'oq 1 

\ 


( s’o^CC) 

c^nc) 


Pi Oi 1 



tpfiC) 

c^no 

= 

p 2 O2 1 



2 ’ 2 ”(C) 

C^k-iiC) 


Pk-1 

Ok-1 1 


^k-iiC) 

C^kiC) J 


\ 

Pk Ok) 


\ 2 ’fe”(C) 


The system at the zeros ( = Q becomes 


'ao - Cj 1 

Pi oi - Cj 1 

P2 0 L 2 - Cj 1 


\ 


2’o"(0) 

^liQ) 


= 0 , 


0 <j<k. 




Pk-l 


Oik-1 Cj 

Pk 


1 

Ok ~ Cj I 


^k-iiQ) 

\^kiCj)J 


( 22 ) 


( 23 ) 


( 24 ) 


Hence, {Cjjj^Q are eigenvalues of the tridiagonal matrix ([ 19 ]). 


□ 


B. Factorization Property 

Another interesting property of the polynomials {S’fc (C)}fc=o) ii^^t to being an orthogonal sequence, is that when the 
parameter n takes positive i nteg er values the polynomials exhibit a factorization property similar to the Bender-Dunne 
orthogonal polynomials [ij, [l 3 , [3 • 

Clearly, The factorization occurs because the third-term in the recursion relation (IlSp vanishes when fc = n + 1 , so that 
all subsequent polynomials have a common factor (Pifi_^_i(C) called a critical polynomial. To illustrate this factorization 
property, we examine it in the case of n = 1 : 

2’l(C) = C, 

®2 (C) = “ 02,iC + 02.002.2, 

^ 3^0 = (C- 2o2.i - 2 a 3 ,i)!p 2 HC), 

^4(0 = (602,1 “ 302,002,2 + 1802,103,1 + 1203,1 - 602,003,2 - 502 ,iC “ 803,iC + C ^) ^2 (C)- 
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Here, the critical polynomial P^iO factorize every other polynomial -P^+2(C): k = 1,2, .... As another illustration, we 
note for n = 2 : 

= c, 

(C) = “ 02,iC + 202,002,2, 

^3 (C) = “ 402 , 002 , 102,2 “ 402 , 002 , 203,1 + (202,1 402 , 002,2 + 2o2,103,1 + 202,003,2) C “ ( 302,1 + 2 o 3 ,i)C^ + C^, 

^pKC) = (C-302,1 -603,i)2’|(C), 

^5 (C) = (1202,i - 4o2,0O2,2 + 6O02,103,1 + 72 o 3 1 - 12o2,003,2 - ( 702,1 + 18 o 3 ,i)C + 

with (C) being a critical polynomial factorize every other polynomial Pk+siOt^ =1,2 ,.... Indeed, all the polyno¬ 
mials ^k+n+iiO, beyond some critical polynomial !P"+i(C), factored into the product 

!pfc”,„.i(C) = Qfe(CX.i(C), fc = o,i,.... ( 25 ) 

Interestingly, the quotient polynomials {Q/l(C)}fc>o form an infinite sequence of orthogonal polynomials. To prove 
this claim, we substitute (j 25 p into (TT^ where we re-index the polynomials to eliminate the common factor 21 ^+i(C) 
from both sides. The recurrence relation (jlSI) then reduces to a three-term recurrence relation for the polynomials 
{ 2 fc(C)}fc>o that reads 

2 fc(C) = iC - (k + n){{k + n - l)a 3 ^i + 02,1)) Qk-iiO - (fc - l)(/c + n)a2,o((A: + n - 1)03,2 + 02,2)2^2(0, ^ ^ 1 , ( 26 ) 

initiated with Q"i(C) = 0, and 2 o(C) = 1 - Hence, the quotient polynomials Q ^{0 ^ 1 ®° form a new sequence of 
orthogonal polynomials for each value of n. The first few polynomials of Q^iO can to explicitly evaluate to reads 


2i(C)-C“(^+ 1)(«2 ,i + ^^a3,l), 

22 (C) = “ (203,1(1 + "o)^ + 02,1 ( 2 n + 3 ))C + {n+ 2 ) [(1 + 'o)(a2,i + 03,in)(02,1 + 03,1-1- 03,in) - 02,0(02,2 + 03,2 + 03,20)]. 


The Christoffel-Darboux formula for this sequence of orthogonal polynomials reads 

^ Q^iOQ-iC) 2 Li(C) 2 ^(C')- 2 ^(C) 2 Li(C') 

j=o jl(n+ 2)^(02,003,2)'’ (n + 1 + k\(n+ 2)^(02,003,2)'" (n + 1 + (C “ C') 

and as ^ C, 

f _ ( 2 "(C))' _ [QUiOYQUO-miOYQUiO 

j=o j!(n + 2)^(02,003,2)'’ (n + 1 + k\{n + 2)fc(o2,003,2)'" ^n + 1 + ^77)^ 


( 27 ) 


( 28 ) 


C. Moments and weight functions 


From the general theory of orthogonal polynomials, there exists a certain weight function 'JT(^) normalized as 




( 29 ) 


for which 


f ¥i:iO^Y'i0^iOdC = PkPk'Skk', k,k' = 0 ,l,...,n, ( 30 ) 

Js 

where S is the support of the measure ‘W’(C) d( on the real line and Skk' is Kronecker’s symbol. 

Theorem III. 3 . For the monic orthogonal polynomials { 5 ’^((')}fe=o 

^{OdC = plSkj, 


Vj < k. 


( 31 ) 
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Proof. By writing p^,(C) = ^j=o , where ak'-k' = 1 , equation (l 22 l) yields 


hj ^ 

Y^OLk'-j / C^^k{ 0 '^iOdC = PkPk'Skk', k,k' = 0 , 1 ,... ,n, ak';k' = i- 
j=o -'•5 


( 32 ) 


For fc' = 0 , 1 ,..., n, equation ( 1 ^ generate a linear system that reduce to 




/ 1 

0 

0 

0 . 

• O') 

-1 


J^iP^iCWiOdC \ 


'+l;0 

1 

0 

0 . 

. 0 


' Podko ^ 

fsC^kiCWiOdC 


'+2;0 

«2;l 

1 

0 . 

. 0 



fsC^^nCWiOdC 

=/ 

«3;0 

'+3;1 

“3; 2 

1 . 

. 0 

X 

P2dk2 

f^etpjficwiOdc 


<■,0 

aj;i 

<+4.2 

<+4.3 • 

. 0 


Psdks 

j.e'tp^icwiOdCj 


K^k’-O 



<';3 • 

• h 

(fc'+l)x(/i;'+l) 

... 

Pin 


The inverse of the lower triangular matrix is again lower triangular matrix with ones on the main diagonal, 
argument concludes our assertion (HB. . 


( 33 ) 


This 

□ 


From Theorem IIIL 31 we also conclude that 

^^(c)dc = i, f^^j:{onod(: = o ik>i), f^^^{o^j:{0'n^iOdc = o, io<j<k<n). (34) 

Theorem III. 4. The norms of all polynomials with k>n+l vanish. 

Proof. The norms pk can be evaluated using the recurrence relations (fT^ after multiplying throughout by C^~^' 4 t’(C) 
and taking integral over we obtain, because of Theorem IIIL 31 a simple two-term recursion relation for the squared 
norm 

Gk = fc(A:-l-n)a2,o((fc-1)03.2 + 02,2)0^-1, 

with a solution given by 


Gk = | 9 feP = k\ (O 2 ,o 03 , 2 )''(-n)fc 1 

^ , 
la 3 . 2 /fc 


( 35 ) 

which is equal to zero for A: = n+ l,n+2,.... 



□ 

As a result of Theorem IIIL 41 we note 




f [^kiOf no dC = k\ (a2.oa3.2)"(-n)fe (, 

Js Va3.2/fc 

0 < k < n, 


( 36 ) 

/'C 2 ’fe”(C) 2 ’fe-i(C)‘W^(C)dC = fc!(a 2 .oa 3 , 2 )"(-n )4 —] . 
Js \a 3 , 2 )k 

k=l, 2 ,..., 


( 37 ) 

f C [TkiOT no dC = k({k- 1)03.1 -1- 02.i)fc! (02.003.2 
J s 

:)'=(-) , 
\a 3 . 2 ;fc 

fc = 0 ,1,2 ,.... 

( 38 ) 


Although we don’t have a general formula for evaluating the moments 

Pk= [ C'^WiOdC, k = 0 ,l,...,n, 
Js 


(39) 













it is possible to compute /ife recursively using Theorem IIII.3I and the recurrence relation (jl 6 |) . It is not difficult to 
show that 


f lT(C)(iC=l, f CW^(C)<^C = 0, f CW{Qd( = pI = - na 2 , oa 2 , 2 , 02,002,2 < 0, 


and for all fc > 2 , 


For example, 


^Cfffc(C)w^(C)dC = o. 

J^C^W{C)dC = a 2 ,ipl, J^C*W{C)dC=pl + (a 2 p-na 2 ,oa 2 , 2 )pl, and so forth. 


(40) 


(41) 


(42) 


Theorem III. 5. The norms of all polynomials reads 


- fc! (n + 2)k ( 02 , 003 , 2 )^ j n + 1 + 


( 1 “2,2'l 

I n + 1 +- I 

V “ 3 , 2 /j 


(43) 


Proof. The proof follows by multiplying the recurrence relation (|26l) by ^W^(C) integrate over (. This procedure 
implies the two-term recurrence relation 


Qf = k{k + n+ l)a 2 ,o((fc + 0 ) 03,2 + 02 , 2 ) 0 ^ 1 , ^ > 1 

where Qf = fg \ Qff.{C,)\^W{z)dz = fg C’^Q^{Qdf. with a solution given (|43)) . 


(44) 

□ 


IV. EXAMPLE: HYPERBOLIC POTENTIALS 

Recently, Downing [I^ studied a class of potentials characterized by two physical parameters Vq and d shaping the 
potential depth and width respectively 


, sinhy/d) 

V '' U 6 . 


-00 < X < CO. 


cosh°(x/d) 

In this section, we consider first a more general class of potentials given by [if 


(45) 


V{x) = -Po 


sinh^(a:/d) 
cosh‘s (x/d) ’ 


-00 < X < 00 , Vo > 0 , p> 0 , q > p. 


(46) 


The condition p > 0 required to avoid a singularity at x = 0 and he condition p < 5 is to ensure the existence of 
a bound-below potential support the presence of at least one bound state. Indeed, it is well known M- Given a 
non-positive potential P < 0 with 


f V{x)dx < 0 , 

J — 00 


(47) 


then there exists a bound state (independent of the potential depth and range) with energy E < 0 for the Hamiltonian 
H = + V(x). For the potential (H51) . it is not difficult to prove that 


/■ 




p>0, p<q. 


The stationary one-dimensional Schrddinger equation for a nonrelativistic particle of mass m and energy E in a 
hyperbolic potential V (x) written as follows 


— —h( sinhP(x/d) 

2m dx'^ ^ °cosh'^(x/d) 


4>{x) = Eip{x), 


Ro>0, 


-00 < X < 00 . 


(48) 
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can be written using the substitution z = x/d as 

= sd^-fpiz), -oo<z<oo. ( 49 ) 

ax^ cosn^( 2 ) 

where Uq = 2mVo/ft^ and e = 2mElh‘^. Upon making the change of variable 77 = l/cosh^(z) such that the domain 
-00 < z < 00 maps to 0 < ?7 < 1 , equation (H51) yields 

4 ? 7^(1 - 77 ) ^ _ 0^21 d'tpiri) ^ tpiv) = 0; 0 < 77 < 1 , p > 0, q>P- (50) 

dr]^ drj 

It is straightforward to verify that 77 = 0 is a regular singular point with exponents {±d^^/2} where £ < 0 and the 
singular point 77 = 1 is also a regular singular point with exponents s = { 0 , 1 / 2 } utilizing the symmetric and anti¬ 
symmetric solutions [H, respectively. From the classical theory of ordinary differential equation, we know that if 
77 = 00 is a regular singular point as well, then the differential equation (I50|) is a hypergeometric-type equation with 
exact solutions expressed in terms of the Gauss hypergeometric function 2 ^ 1 - This claim can be be easily verified in 
the case of p = 0 and q = 2 where the differential equation (l50l) turn to be Fuchsian with the three regular singular 
points at 77 = 0 , 1 , and 00 . 

The point 77 = 00 is an irregular singular point only if p > 2 and q> p. A reason for V{x) = -Vq sinh^(? 7 )/cosh®( 77 ) to 
be a quasi-exactly solvable potential that will be our focus because of the recent interest. In this case, equation (EOl) 
reads 


77^ (1 - 77) + rj 


2 . 


-H Uo (1 - ’nf) '0(d) = 0- 


Because of the asymptotic behavior at the singular points, the exact solutions takes the form 
7 /( 77 ) = ?7^/^(I-77)®e““''/^/(77), /3, a>0, s = 0,1/2, 0<77<1, 

we obtain, for the function /{rf), the following differential equation 


(51) 


(52) 


/"(d)+ 


1 -H 4s I+ 13 

■ +-a 


2(d-l) d 


,/'(d) + 


-a + (3 + (3^ + As - Aas + Aj3s ^ -2a - (3 - 2af3 - - As - 4/3s + a^ 


4 ( 77 - 1 ) 


477 


/(d) = 0 
(53) 


where s = 0,l/2, /3 = d\f^ and a = d\/Uo. Equation (1531) can be written as 

4^(77 - + [-Aap^ -t- (6 + 4a + 4/3 + 8s)p - 4(1 + /3)] /'(p) 

+ [(a(a - 2/3 - 4s - 3)p - + 2a(l + P) + {1 + /3)(/3 + 4s)] /(p) = 0 (54) 

Using Theorem lII.il the necessary condition for the differential equation (I54|) to admit the polynomial solutions is 

a - 2/3 = 4n + 4s + 3 


that implies the following eigenvalue spectra: 




(477 + 4s H- 3 - d\/Uo) 
4d2 


d^ Uq > (477 + As + 3) ^. 


Again using Theorem lII.il the sufficient condition follows by the determinant 


Po 


0 

0 

0 ... 

0 

0 

0 

7i 

Ml 

^2 

0 

0 ... 

0 

0 

0 

0 

72 

M2 

^3 

0 ... 

0 

0 

0 

0 

0 

73 

M3 

<54 ... 

0 

0 

0 

0 

0 

0 

0 

0 ... 

7n-l 

l^n-1 

dn 

0 

0 

0 

0 

0 ... 

0 

In 

l^n 


(55) 


(56) 
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where 

Mi = + 2a(2j + /3+ 1) + (2j +/3 + l)(2j +/3 + 4s), 7 ^- = 4a(l - j + n), 6j = -4j'(/3 +j)- 

The exact solutions are given explicitly by the equation 

t (-if. 


(57) 


where the coefficients with ( = np = (1 + d^J-Sn) (d^-Sn + 4s) + 2d{l + d-/^) \/Uo-d‘^Uo are evaluated using 

the following three-term recurrence relation 

2’fe”i(C) = (C-2A:(l + 2a+2/3 + 2fc + 4s))fP,"(C)-16fc(fc-n-l)a(/3 + A:)2>,%(C), 2’”i(C) = 0, 2’o”(C) = 1, (58) 

with the square norm given by 

/'[!Pfc"(C)]" 4M(C)dC=2""fc!a"(-n)fc (1 + /3),, 0 < fc < n. (59) 

■J s 

The first few solutions are: 

• For n = 0, we have 


(4s + 3 - (i\/t/o) 

^0 =- -^2 -> foiv) = 1 , 

subject to the root finding of the equation 

7d2[/o-12d^Cio-3(4s+l)(4s+5) = 0 . where Uq > {As + 5)"^jd^. 
The critical polynomial reads 

(a,/?; s) = -a^ + 2a(l + /?) + (1 - 1 - /3)(/3 + 4s). 

For n = 1, we have 

(4s + 7- (^^/C7o) .f_ 2 ^ |q;^ - 2a(l -t- /3) - (1 4 - /3)(/3 + 4s) 


(60) 

(61) 


£1 = 


4d2 


fiiv) 


4(1+ /3) 




(62) 


(63) 


subject to the root finding of the equation 

- 4(2 -I- /3)a^ + 2(3 + 5)3 + 13'^- 4(2 + I3)s)a^ + 4(1 + /3)(3 -i- 4/3 + + 4n - 1 - 4(3 + /3)s)a 

+ (l-i-/3)(3 + /3)(/3 + 4s)(2 + /3 + 4s) = 0 

where a = d\/Uo, (3 = (7 + 4s - d\/Uo)l2 and d^Uo> (7 + 4s)The critical polynomial in this case reads 

(P 2 {^T13] s) = 16q;(1 + /3) + {—(x^ + 2o;(l 4- /3) 4- (1 4- 13) {)3 + 4s)) (— 4- 2o:(l 4-/3)4-(l + /3)(/3 4' 4s) ~2(34-2o;4-2/34- 4s)). 

For n = 2, we have 

(4s 4-11 — d\/Uo) 


£2 = 


/2(m) = 1- 


4d2 

c? - 2a(l 4- /3) - (1 4- /3)(/3 4- 4s) 
4(1+ /3) ’ 


32q;(1 4- /3) - (6 4- 4- 2a(l - /3) 4- 4s 4- /3(3 - /3 - 4s))(-a^ 4- 2a{l 4- /3) 4- (1 4- /3)(/3 4- 4s)) , a\ 

"-’ * ’ 

subject to the the root finding of the equation 

32o:(2 4 - /3)(— 4 - 2o:(l 4 - /3) 4 - (1 4 - /3)(/3 4 - 4s)) 

4- {—cP' + 2(x{5 + p) 4- (54-/3)(44-/3 + 4s))(32o:(l4-/3) 4- {—cP + 2o;(l 4- /3) 4- (1 4- /3)(/3 4- 4s) ) 

X {-a^ + 2a(3 4- /3) 4- (3 4- /3)(2 4- /3 4- 4s))) = 0 (65) 

where a = d\/Uo, /3 = (11 4- 4s - d\/Uo)l2 and d'^Uo > (11 4- 4s)^. The critical polynomial reads in this case as 

2>3 (a, /3; s) = 32a(l 4- /3)(-a^ 4- 2a(l 4- /3) 4- (1 4- /3)(/3 4- 4s) - 4(5 + 2a + 2P + 4s)) 

4 - (2q:(1 4 - P) — cP 4 - (1 + P){P + 4s))(32o;(24- P) 4 - {—cP + 2q:(1 4 - P) + (1 4 - P){P 4 - 4s) — 2(3 4 - 2 cr 4 - 2/3 4 - 4s)) 

X {-a^ + 2a(l 4- /3) 4- (1 4- /3)(/3 4- 4s) - 4(5 4- 2a 4- 2/3 4- 4s))) 
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V. CONCLUSIONS 

We have given the sufficient and necessary conditions that ensure the existence of the polynomial solutions of a 
general confluent Heun equation. It is indeed exciting to know that these polynomial solutions are the source of 
finite sequences of orthogonal polynomials. The number of the orthogonal polynomials in each sequence equal to the 
degree of the polynomial solution of the confluent Heun equation. Although we discussed several properties of these 
orthogonal polynomials, a general theory that explain the zeros and evaluation of the discrete weight function will 
be of great interest. We hope the present study will stimulate further research on the theory of finite orthogonal 
polynomials. 
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